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ABSTRACT 

We present new optical, near-IR, and mid-IR observations of the young eruptive variable star 
V1647 Orionis that went into outburst in late 2004 for approximately two years. Our observations 
were taken one year after the star had faded to its pre-outburst optical brightness and show that 
V1647 Ori is still actively accreting circumstellar material. We compare and contrast these data 
with existing observations of the source from both pre-outburst and outburst phases. From near- 
IR spectroscopy we identify photospheric absorption features for the first time that allow us to 
constrain the classification of the young star itself. Our best fit spectral type is M0±2 sub-classes 
with a visual extinction of 19±2 magnitudes and a K-band veiling of ik ^1.5±0.2. We estimate 
that V1647 Ori has a quiescent bolometric luminosity of ^^9.5 Lq and a mass accretion rate of 
^lxlO~^ Mq yr^^. Our derived mass and age, from comparison with evolutionary models, are 
0.8±0.2 Mq and < 0.5 Myrs, respectively. The presence towards the star of shock excited optical 
[S II] and [Fe II] emission as well as near-IR H2 and [Fe II] emission perhaps suggests that a new 
Herbig-Haro flow is becoming visible close to the star. 

Subject headings: stars: individual (V1647 Ori) - Reflection nebulae - Accretion, accretion disks 



1. INTRODUCTION 

A significant event in the recent history of star 
formation studies occurred in late 2003 when the 
young pre-main sequence star V1647 Orionis went 
into outburst. This eruption, and the associated 
one hundred fold increase in optical brightness, re- 
sulted in the appearance of a monopolar reflection 
nebula now known as " McNeil's Nebula" after the 
amateur astronomer. Jay McNeil, who discovered 
it (McNeil 2004). The star and nebula remained 



bright for approximately 18 months before fading 
rapidly over a six month period. V1647 Ori was 
found to be again close to its pre-outburst (Sloan 
Digital Sky Survey, SDSS) optical brightness in 
early 2006. A widely accepted interpretation of 
such eruptions is one involving a heated disk and 
a 'mass dumping' episode from a circumstellar disk 
onto a stellar photosphere. The five magnitude in- 
crease in optical brightness seen in V1647 Ori has 
been attributed to both the addition of significant 
accretion luminosity and a partial dust clearing 
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due to the high- velocity wind (up to 600 km s ^) 
emanating from the star /disk engine (Reipurth & 
Aspin 2004; McGehee et al. 2004). 

Both star and nebula have been widely studied 
from the period soon after the outburst occurred 
through to its fading to its pre-outburst bright- 
ness. In our discussions below, we refer to the 
period before November 2004 as the pre-outburst 
phase (Briceno et al. 2004), the period from 
November 2004 to February 2006 as the outburst 
phase, and the period from February 2006 onwards 
as the quiescent phase of V1647 Ori. Observations 
spanning all spectral regimes, from x-ray into the 
optical to the near-IR and on to the far-IR/sub- 
mm/radio, have been published by many authors 
in a large number of papers. Readers are referred 
to these for detailed background information on 
V1647 Ori and McNeil's Nebula. These papers 
are, in alphabetical order: Abraham ot al. (2004 - 
the pre-outburst infrared characteristics), Acosta- 
Pulido et al. (2007 - optical and near-IR imaging 
photometry, near-IR spectroscopy and polarime- 
try during the outburst, henceforth AP07), An- 
drews, Rothberg, & Simon (2004 - mid-IR and 
sub- mm observations during the outburst), As- 
pin et al. (2006 - an historical study of previ- 
ous outbursts), Briceno et al. (2004 - the opti- 
cal outburst history, pre-outburst and during the 
outburst), Fedele et al. (2007a,b - optical and 
near-IR properties during the outburst), Gibb et 
al. (2006 - near-IR spectroscopy during the out- 
burst), Grosso et al. (2005 - x-ray observations 
during the outburst), Kastner et al. (2004 - x-ray 
emission pre-outburst and during the outburst), 
Kastner et al. (2006 - x-ray evolution during the 
outburst), Kospal et al. (2005 - optical photome- 
try during the outburst), McGchcc et al. (2004 - 
optical and near-IR photometry pre-outburst and 
during the outburst), Mosoni et al. (2005 - mid-IR 
intcrfcromctric observations during the outburst), 
MuzeroUe et al. (2005 - Spitzer observations dur- 
ing the outburst), Ojha et al. (2004 - near-IR 
imaging during the outburst), Ojha ct al. (2006 
optical photometry and spectroscopy, and near-IR 
imaging during the outburst), Reipurth & Aspin 
(2004 - optical and near-IR imaging photometry 
and spectroscopy during the outburst), Rettig et 
al. (2005 - high spectral resolution near-IR ob- 
servations during the outburst), Scmkov (2004 - 
optical photometry during the outburst), Semkov 



(2006 - optical light curve during the outburst), 
Tsukagoshi et al. (2005 - mm continuum obser- 
vations during the outburst), Vacca, Gushing, & 
Simon (2004 - near-IR spectroscopy during the 
outburst), Vig et al. (2006 - radio observations 
during the outburst), and Walter et al. (2004 - 
optical photometry and spectroscopy during the 
outburst). 

The significance of eruptive outbursts, similar 
to the one undergone by V1647 Ori lies in the fact 
that these periods are considered times when mass 
accretion rates increase dramatically (Hartmann 
& Kenyon 1996, henceforth IIK96). During peri- 
ods of outburst, mass accretion rates have been di- 
rectly observed to increase by up to several orders 
of magnitude from a pseudo-steady-state value of 
~10-9 *o -7 Mq yr-i for classical T Tauri stars 
(Hartmann et al. 1998) to ~10-^ *° M© yr"! 
for EXor variables, and as high as ^lO^** Mq yr^^ 
for the extremely energetic FUor variables (Hart- 
mann fc Kenyon 1985). FUors and EXors are so 
named after the prototypes of their classes, namely 
FU Orionis and EX Lupi, respectively (Herbig 
1989). 

Whether all young stars undergo EXor and/or 
FUor eruptions and, if so, whether these two 
types of outbursts events are fundamentally the 
same (on perhaps different scales) or difi'erent (in 
terms of the mechanism involved and the trigger- 
ing of the eruption) is still a matter of debate. 
EXor eruptions are distinct from FUor eruptions 
in that they are much shorter lived, months to 
years rather than decades to centuries, and have 
been empirically determined to be repetitive (Her- 
big 1989; 2007). For example, EX Lupi itself has 
gone into outburst on several occasions, with the 
last being in 1993 94 (Lchmann et al. 1995, Her- 
big ct al. 2001). V1647 Ori is likely another ex- 
ample of an EXor since its recent activity lasted 
only two years and it was previously observed in 
outburst from 1966-67 (Aspin et al. 2006) on 
archival photographic plates. Mechanisms that 
have been proposed to initiate an eruptive event, 
be it EXor-like or FUor-like, include: thermal disk 
instabilities resulting in a runaway accretion con- 
dition (Bell & Lin 1994), periodic overloading of 
the inner regions of the circumstellar disk and sub- 
sequent magnetic collapse (HK96), and the close 
approach of a companion in an eccentric orbit dis- 
turbing the stability of the inner regions of the 
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disk (Bonnell & Bastien 1992). 

In this paper, we present new optical, near-IR, 
and mid-IR observations of V1647 Ori and Mc- 
Neil's Nebula taken approximately one year after 
the source reached its pre-outburst optical bright- 
ness. Below, we compare and contrast our obser- 
vations with those taken during the pre-outburst 
and outburst phases of the eruption and attempt 
to better understand the nature and characteris- 
tics of the underlying young star and its circum- 
stellar environment. 

2. OBSERVATIONS & DATA REDUC- 
TION 

In order to study the on-set of the quiescent 
phase of V1647 Ori, we undertook an observing 
campaign involving optical, near-IR, and mid-IR 
instrumentation. Table 1 relates the log of obser- 
vations obtained together with associated photo- 
metric values derived from those datasets. All the 
optical and mid-IR observations were acquired on 
the "Frederick C. Gillett" Gemini North 8- meter 
telescope on Mauna Kea, Hawaii during the period 
UT February 13 to April 5 2007. At this time, 
VI 647 Ori had returned to its pre-outburst quies- 
cent brightness for about one year. Additionally, 
near-IR images and spectroscopy were acquired in 
February 2007 using the University of Hawaii 2.2m 
and the NASA IRTF 3.5m telescopes, respectively, 
both also located on Mauna Kea. 

2.1. Optical Imaging & Spectroscopy 

Optical imaging data using Gemini North and 
GMOS-N (Davies et al. 1997) were obtained un- 
der Gemini queue program ID GN-2007A-Q-33 on 
UT February 22, 2007. Images in the SDSS g', r', 
i', and z' filters were acquired in queue mode with 
an exposure time of 10 minutes per filter. The 
data were reduced with the Gemini IRAF data re- 
duction gmos package (vl.9) using standard proce- 
dures for bias subtraction, flat-fielding (using mas- 
ter twilight fiats), cosmic ray rejection, and image 
combination. Photometric calibration was per- 
formed using a photometric sequence of six field 
stars from SDSS images of the region. Since this 
is a region of active star formation with the likeli- 
hood of some stars being intrinsically variable, we 
considered using as many as six field stars as cal- 
ibrators to be an acceptable way to ensure insen- 



sitivity to such temporal changes. A detailed dis- 
cussion of the calibration sequence definition and 
properties will be presented in a subsequent pa- 
per on our two year outburst phase Gemini/IRTF 
optical, near-IR, and mid-IR monitoring program 
on V1647 Ori (Aspin & Reipurth, in prep.). The 
GMOS CCDs were binned 2 x 2 to match the natu- 
ral seeing of ~0'.'55. A 'true-color' image of the re- 
gion containing V1647 Ori, made from these data, 
is shown in Fig. [TJ Also shown is a similar image 
taken some three years earlier during the outburst 
phase (February 2004). 

Additionally, optical spectroscopy of V1647 Ori 
was obtained on UT February 21, 2007, again 
using Gemini North and the GMOS-N in queue 
mode. We utilized the R400 grating with Ac=7000 A 
and a 0'.'75 wide long-slit. This resulted in a spec- 
tral resolution of R'^1280 and wavelength cover- 
age of ~5000-9150 A. Seven 1800s exposures were 
taken giving a total on-source exposure time of 
3.5 hours. V1647 Ori was moved along the slit 
between each observation and, at the time of our 
optical spectroscopy, it had an r' magnitude of 
23.26±0.15. The data were again reduced using 
the Gemini IRAF gmos package (vl.9) using flat- 
field and arc images from the facility calibration 
unit, GOAL. In order to obtain the best sky-line 
subtraction in the far-red, we subtracted pairs 
of reduced 1800s images before residual sky-line 
subtraction and optimal point-source extraction. 
These steps were performed using the Starlink 
FIGARO package. The seven extracted spectra 
of V1647 Ori were inter-compared to eliminate 
spurious noise spikes before combination. The 
combined spectrum was flux calibrated using ob- 
servations of the spectrophotonietric standard star 
Feige 34. Fig. [2] shows the final optical spectrum 
of VI 647 Ori with all significant atomic lines iden- 
tified. 

2.2. Near-IR Imaging & Spectroscopy 

We obtained near-IR images of V1647 Ori on 
UT February 13, 2007 on the University of Hawaii 
2.2 m telescope on Mauna Kea using ULBCAM 
(Hall et al. 2004). This instrument consists of 
four Hawaii-2 1-5 fj,m arrays but operates in the 
J and H band only due to thermal background is- 
sues. The field of view of the camera is 17'xl7' 
with a pixel scale of 0'.'25. A standard "Mauna 
Kea" H-band filter was used for the observations 
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and the total on-source exposure time was 6 min- 
utes. Target images were dithered and median 
filtered to obtain sky images and flat-field images 
were taken by observing the inside of the telescope 
dome with lights on and off. Data reduction was 
performed using the Gemini IRAF quire data re- 
duction package (vl.9). For photometric calibra- 
tion we used 2MASS sources within the field ob- 
served. Fig. [3] shows the region immediately sur- 
rounding V1647 Ori from these H-band images. 
At the time of the observations, V1647 Ori had 
an H-band magnitude of 11.96±0.1 and is labeled 
in Fig. [3] together with LkHa 301 and HH 22-IRS 
(MuzeroUe et al. 2005). 

We used the NASA IRTF and the near-IR spec- 
trometer SpeX (Rayner et al. 2003) to obtain 
near-IR 1-4.1 /im spectroscopy of V1647 Ori on 
UT February 22, 2007. We observed using the 
SXD and LXD grisms with total exposure times 
of 1600 and 1200 seconds, respectively. The data 
were reduced using the SpeXtools software pack- 
age (Gushing, Vacca, & Rayner 2004). An ob- 
servation of an AO V star close to the airmass of 
the target was used as a telluric calibrator and 
applied using the template fitting method defined 
in Vacca, Gushing, & Rayner (2003). V1647 Ori 
was nodded along the O'.'S slit in an ABBA pat- 
tern to allow good sky definition and subtrac- 
tion. The final reduced and combined spectrum 
was flux calibrated using the H-band magnitude 
(H=11.96) from our ULBGAM imaging taken 9 
days prior to eliminate slit-losses during the vari- 
able seeing conditions encountered. The conver- 
sion of H magnitude to fiux was performed us- 
ing the Tokunaga & Vacca (2005) values, for 
which a zero H magnitude star has a fiux of 
Fa=1.18 10"^ W m-2 ^m-i. The combined SXD 
and LXD spectrum of V1647 Ori is presented 
in Fig. m where all major absorption and emis- 
sion lines and bands are identified. An expanded 
view of the 2.1-2.35 /im region of the spectrum of 
V1647 Ori is shown in Fig. [5] 

2.3. Mid-IR Imaging & Spectroscopy 

On UT March 18, 2007 we obtained Gemini 
North Michelle images of V1647 Ori in queue mode 
using the N' (A iim) and Qa (A c=18.5 /xm) 

filters. The total on-source exposure time was 
100 seconds in each filter. The chopped and nod- 
ded images were processed using the Gemini IRAF 



midir package (vl.9), specifically task msreduce. 
Photometric calibration was performed using ob- 
servations of the bright photometric standards 
56 Ori and HD37160. Since V1647 Ori appeared 
in both images as a point-source with N' and Qa 
band full width half maximum seeing of 0'.'48 and 
0'.'66, respectively, they are not explicitly shown 
here. From these images, however, we find that 
V1647 Ori had flux values of 0.23 Jy at 11.2 /xm 
and 0.44 Jy at 18.5 /im. 

Additionally, mid-IR spectroscopy covering the 
8-14 /im wavelength regime (the N-band) was ob- 
tained on Gemini North using Michelle with the 
'low N' grating on UT April 5, 2007. The slit 
used was 0'.'4 wide and resulted in a spectral res- 
olution of R^200. The on-source exposure time 
was 400 seconds and chopping along the long-slit 
(±10") was employed. A spectroscopic standard 
star (BS2473, spectral type G8 lb) was observed 
immediately after the science observation and used 
to remove telluric features from the target spec- 
trum. Spectroscopic reduction was performed us- 
ing the Gemini IRAF midir package (vl.9) to give 
a reduced 2D spectral image, and the Starlink FI- 
GARO package was used for optimal spectroscopic 
extraction and telluric correction. 

3. RESULTS 

3.1. The Optical Appearance 

In Fig. [1] we compare V1647 Ori and McNeil's 
Nebula on images taken in February 2004 and 
three years later in February 2007. The images 
are derived from three optical passbands, the left 
image (February 2004) uses the g', r', and i' im- 
ages while the right image (February 2007) uses 
the r', i', and z' data. Note that the optical r'- 
band fiux from V1647 Ori is ^170 x smaller in 
February 2007 than in February 2004, the latter 
being 2-3 months after the initial outburst. We 
also note that in February 2007, some 12 months 
after V1647 Ori had returned to its pre-outburst 
brightness, we barely detect the refiection nebu- 
losity and now see the HH 22 region (blue, to the 
north and slightly east of V1647 Ori) much more 
distinctly. Even though the reflection nebula is 
very faint, all the physical features of it are still de- 
tected in the February 2007 image, specifically the 
fan-shaped region extending from the star to the 
north-west, the dark bands crossing the nebula on 
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the eastern side of the monopolar lobe, the right- 
angled feature directly north of V1647 Ori, and 
the small nebulous blob a few arcseconds north- 
east. Since we are observing the nebula in re- 
flection (for example, see the polarization map of 
AP07), the brightness of the nebula should ap- 
proximately mirror the brightness of V1647 Ori 
itself with appropriate delays for light-travel time 
(AP07). We conclude, therefore, that the under- 
lying structure observed in the reflection nebula is 
quasi-static and independent of the eruption and 
brightness changes that have taken place. 

Qualitatively, the scenario for the rapid bright- 
ening and fading of McNeil's Nebula is as follows. 
The initial explosive eruption blew away and/or 
sublimated considerable obscuring dust along both 
our direct line-of-sight and the line-of-sight from 
VI 647 Ori into the nebula (Reipurth & Aspin 
2004; Walter et al. 2004). This was the direct 
result of both the initiation of an intense stel- 
lar wind and the enhancement to the UV/optical 
flux from the star due to the increased accretion 
rate. We note that the bolometric luminosity, 
Lfjoi, of V1647 Ori increased from Lf,,,; ~3.5 L0 
pre-outburst to Li,oi ~50 Lq during outburst (An- 
drews, Rothberg, & Simon 2004) due to the ad- 
dition of enhanced accretion luminosity (Lace)- 
The addition of the Lqcc component and the blow 
out/sublimation of surrounding dust resulted in 
the brightening of V1647 Ori and hence, via re- 
flection, McNeil's Nebula. The reverse is also 
likely true, i.e. the decline of the Lqcc compo- 
nent and the slowing of the intense stellar wind 
resulted in the condensation of dust restricting the 
direct flux and the flux propagating into the neb- 
ula. Dust sublimation was recently discussed by 
AP07 who commented on its reversibility i.e. sub- 
limation during the outburst, and re-condensation 
after fading, and how such behavior was consistent 
with the documented evidence of at least one pre- 
vious outburst (Aspin et al. 2006). 

It has been proposed that the existence and 
shape of McNeil's Nebula is the result of either 
a remnant wide-angle molecular outflow creating 
a cavity in the ambient molecular cloud material 
near V1647 Ori (Reipurth & Aspin 2004), or that 
such nebulae are the result of space motion of the 
young star through the ambient molecular cloud 
(Goodrich 1987). In the former case, an oppos- 
ing outflow lobe would extend into the molecular 



cloud (since the visible lobe is tilted towards us - 
AP07) and hence be hidden from view by cloud 
material. As yet, there is no observational evi- 
dence for an active or remnant molecular outflow 
emanating from V1647 Ori (Lis, Menten, & Zylka 
1999, Andrews, Rothberg, & Simon 2004). Simi- 
larly, there has been no observational determina- 
tion of space motion associated with V1647 Ori. 

3.2. The Optical Spectrum 

In Fig. [2] we present the optical spectrum of 
V1647 Ori taken in February 2007. The top of 
the three panels shows the complete optical spec- 
trum from 5000 A to 9150 A. We note that at the 
time this spectrum was taken, V1647 Ori had an 
optical brightness of r'=23.3±0.1. The spectrum 
seen is that of an object with a red continuum 
super-imposed with a variety of emission features. 
The brightest emission lines are those from the 
Ca II triplet at 8498, 8542, and 8662 A. A pro- 
nounced Ha emission feature, together with weak 
high-order H I Paschen series PIO, Pll and P12 
lines, are also present. Weaker emission lines are 
present also, specifically [S II] at 6716 and 6731 A, 
O I at 5579 and 8449 A, [Fe II] at 8617 A, very 
weak Fe I lines at 7912, 8467, and 8689 A, and 
weak Fe II at 8600 A. The absorption feature at 
7600 A is the telluric O2-A band. The middle and 
bottom panels of Fig. [2] show expanded views of 
the optical spectrum near Ha (middle) and the 
Ca II triplet (bottom). The Li I absorption line at 
6707 A is indicated, however, it is not clear that 
it is really detected. Table 2 shows the equivalent 
widths and line fluxes for the optical (and near-IR) 
emission features. 

3.2.1. Shock- excited fS II] emission 

An interesting aspect of the optical spectrum 
is the presence of weak [S II] emission lines with 
equivalent widths of -16.3 A (6716 A) and -22.0 A 
(6731 A). Such [S II] emission features were also 
detected by Fedele et al. (2007a) in their January 
2006 observations. They interpreted these as pos- 
sible evidence for the presence of a Herbig-Haro 
(HH) object close to the star. We shall consider 
this interpretation further below. 

We have made a careful study of the CMOS 
two-dimensional spectroscopic images to examine 
whether the [S II] emission is present directly on 
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the location of V1647 Ori or is also seen at other 
spatial positions along the long-slit used. Our 
findings indicate that there is no [S II] emission 
at all at any other location along the 90" long- 
slit and we conclude that the shock-excited [S II] 
emission lines are directly associated with the cir- 
cumstellar region of V1647 Ori and are not due to 
foreground or background shocked nebulosity. 

What evidence is there to support the interpre- 
tation that the [S II] lines originate in a shock- 
excited HH object? First, Eisloffel & Mundt 
(1997) considered that HH 23, located 155" north 
of V1647 Ori and close to the axis of McNeil's 
Nebula, originated from V1647 Ori. Second, as 
well as [S II] , the optical and near-IR spectrum of 
V1647 Ori shows weak [Fe II] emission at 8617 A 
and 1.644 /im, and weak H2 v=l-0 S(l) emission 
at 2.122 /im. Some of these lines are often associ- 
ated with shock-excited HH objects (Gredel 1994, 
1996; Hamann et al. 1994). Third, it is becom- 
ing more widely recognized that there very likely 
exists a relationship between the creation of an 
HH object and the outburst phenomenon (Dopita 
1978; Reipurth 1989). 

If we assume that HH 23 was ejected from 
V1647 Ori in perhaps an earlier outburst, then, 
with an inclination of the ejection axis with re- 
spect to the line-of-sight of ~60 degrees (AP07) 
and assuming a typical HH object space velocity 
of ~200 km s~^, HH 23 would have been ejected 
from V1647 Ori approximately 2000 years ago for 
a source distance of 450 pc. However, it is, as yet, 
not definitively proven that HH 23 originates from 
V1647 Ori. In Fig.[6]we show the region containing 
V1647 Ori, HH 22, and HH 23 imaged in narrow- 
band [S II] using GMOS-S on the Gemini South 
telescope on UT October 10 and 11, 2005. This 
deep image clearly shows shock-excited nebulosity 
between V1647 Ori/McNeil's Nebula and HH 23 
which resembles an intermediate bow shock with a 
trailing extended tail. These features were identi- 
fied as HH23 G by Eisloffel & Mundt (1997). Glose 
inspection of Fig. [6l suggests that it is equally 
conceivable that HH23 G and, in fact, the main 
HH 23 knots (A and B) could have originated from 
HH 22. Within the HH 22 A shock-excited nebu- 
losity lies a deeply embedded young Glass I proto- 
stellar source first detected by Spitzer (MuzeroUe 
et al. 2005). We henceforth refer to this source as 
HH22-IRS and identify its location in Fig. [6] with 



a white cross. MuzeroUe et al. (2005) considered 
that HH22-IRS drives the east-west HH 22 colli- 
mated flow and additionally noted that this flow 
is approximately orthogonal to the HH 23 flow. If 
HH22-IRS drives both the HH 22 and HH 23 flows 
then it is perhaps a young binary system in which 
each component drives a (non-aligned) collimated 
flow. Proper motion and radial velocity studies of 
both flows may clarify their origin. 

If the presence of shock-excited optical and 
near-IR emission lines does indicate the presence 
of a new HH object close to V1647 Ori, when was 
this HH flow created? We have investigated the 
possibility that these shock-excited [S II] lines were 
present before January 2006, and hence could have 
been created in the explosive outburst itself. In 
this case, the [S II] lines would have been present 
in spectra of V1647 Ori taken prior to January 
2006 but were perhaps overwhelmed by the then 
100 X brighter continuum flux. We note that in 
January 2006, V1647 Ori had faded significantly 
from it outburst maximum brightness of r'=17.4 
(Reipurth & Aspin 2004). We have simulated ob- 
serving the same line strength [S II] emission fea- 
tures on such an enhanced continuum. In addi- 
tion to considerably reducing the visibility of the 
lines on the high continuum level, a factor 100 x 
increase in continuum flux implies the addition of 
10 X the photon shot noise to the data. Under 
these conditions the [S II] equivalent widths be- 
come 100 X smaller and the noise becomes 10 x 
larger. We measure the root-mean square (rms) 
noise in a continuum subtracted region of the 
V1647 Ori spectrum close to the [S II] lines at 
3.5xl0~^^ W m~^ /im^^ while the peak signal in 
the [S II] lines is 5x10^20 w m'^ /xm"!. This 
implies a signal-to-noise ratio (S:N) of ^^14. If we 
increase the photon shot noise by a factor x 10 the 
S:N drops to ~1.4 and would mean that the [S II] 
lines of the same strength would likely not be de- 
tected during the outburst phase. 

What we can conclude therefore, is that [S II] 
emission was present in our February 2007 spec- 
trum of V1647 Ori as well as the January 2006 
spectrum of Fedele et al. (2007a). We have shown 
that the emission is isolated to the point-source it- 
self but cannot determine whether these emission 
lines have recently formed or if they were present 
in earlier spectra and swamped by the enhanced 
continuum from V1647 Ori. Again, if we assume 
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an HH space velocity of ~200 km and an in- 
clination of ~60 degrees, a new HH object would 
take ^13 years to move ~1" (for a distance of 
450 pc). 

The fact that [S II] emission lines are present, 
however, implies that shocks exist in the imme- 
diate vicinity of the accreting young star. Since 
the [S II] 6716 A/6731 A hue ratio is only weakly 
dependent on temperature, it can be used as a di- 
agnostic of the electron density in the [S II] emit- 
ting region. Using the IRAF nebular . temden 
routine and our February 2007 observed line ra- 
tio of ~0.8, we estimate that the electron den- 
sity, Tie, was ^^1400 cm""^ for an assumed elec- 
tron temperature, Te, of lO** K (a value consistent 
with mild shock excitation). If we vary Tg over 
the range 5000-20000 K (the temperature range 
over which [S II] emission can occur), a ratio of 
0.8±0.1 gives a range of of -1000-1800 cm'^. 
This derived value is a factor '^14x lower than the 
[S II] critical density of ^^2x10^ cm"'^ (Bacciotti 
2002). Performing the same analysis for the Fedele 
et al. (2007a) data from January 2006, we find a 
[S II] line ratio of ~0.6±0.1 and an electron den- 
sity value of ne~4200 cm~^ with a range for valid 
Te values of 3200-5300 cm'^. The Fedele et al. 
(2007a) value is perhaps a little larger than than 
found in February 2007 although the errors on the 
ratios are significant. A difference in this line ra- 
tio would mean that conditions in the region where 
the shock-excitation occurs have evolved over the 
intervening year and could be interpreted as ei- 
ther physical electron density changes or changes 
in the local ionization fraction. There is some ev- 
idence that in HH flows, higher electron densities 
may occur closer to the exciting sources (Cohen & 
Jones 1987; Biihrke, Mundt, & Ray 1988; Heath- 
cote & Reipurth 1992; Beck et al. 2007). If the 
above ratio decrease is real, then we can speculate 
that perhaps the shocked region has moved further 
away from the star. 

3.2.2. Ha emission 

The equivalent width of the Ha emission in 
our spectrum is Wa=-124±10 A, the large uncer- 
tainty being due to the very weak continuum flux 
present at the wavelength of Ha. The observed 
Wa value is considerably larger than that last re- 
ported: Ojha et al. (2006) found Wa=-40.75 on 
UT September 28, 2005. This is also likely to be 



a result of the weak observed continuum flux. 

The Ha emission line flux in our spectrum is 
Fa=:1.0x 10"^^ W m~^ and we can compare this 
directly to the outburst phase Ha line fluxes pre- 
sented in AP07, Ojha et al. (2006), and Walter et 
al. (2004). AP07 plot their own Ha hne fluxes, 
together with those of Ojha et al. (2006) and 
Walther et al. (2004), against time (their Fig. 12) 
from soon after the outburst was detected in early 
2004 through late 2005. Their plot shows a declin- 
ing Ha line flux during the outburst period, from 
a peak in early 2004 of -2.0x10-^^ W m'^ to 
--4.0x10"^^ W m"2 in late 2005. Our quiescent 
phase value from early 2007 is lower still by a fac- 
tor x4 and follows the outburst trend. As pointed 
out by AP07, H I emission lines are often directly 
related to the accretion process. This trend could 
represent declining accretion over the outburst to 
quiescent phase period although as we shall see 
below, changes in visual extinction could also af- 
fect the observed fluxes. We defer a discussion of 
derived accretion rates during these phases to Sec- 
tion 3.7.3 where we consider the near-IR spectrum 
of V1647 Ori in relation to the models of MuzeroUe 
et al. (1998, 2001). 

We additionally note that the Ha line profile 
appears somewhat asymmetric in shape. Specifi- 
cally, the red wing appears truncated with respect 
to the blue wing and with respect to the profile of 
other emission lines (e.g. Ca II). Such a red-wing 
truncation has not, to our knowledge, been seen 
before in V1647 OriQ We can perhaps attribute 
this asymmetry to redshifted absorption from in- 
falling cool neutral gas although higher spectral 
resolution optical data would be required before a 
conclusive decision can be made. 

3.2.3. Ca II triplet emission 

In the quiescent-phase spectrum shown in 
Fig. H the Ca II triplet hues at 8498, 8542, and 
8662 A have equivalent widths, Wa, of-30.6±l A, 
-35.4±1 A, and -34.2±1 A, respectively. The ra- 
tio of these lines, specifically, ^^1. 7:2:1. 9, is some- 
what different from that found by Walter et al. 
(2004) in spectra from the outburst period taken 
between February and April 2004. They reported 



A recent high spectral resolution Keck II/NIRSPEC near- 
IR spectrum of V1647 Ori taken by us shows a similarly 
truncated Br7 emission line. 
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a Ca II triplet ratio of ^--^2:2:1, a pattern consid- 
ered unique to T Tauri and Herbig Ae/Be stars 
(Herbig & Soderblom 1980; Hamann & Persson 
1990), and inconsistent with the expected ratio 
from (optically thin) transition strengths (gf) i.e. 
1:9:5. The 2:2:1 line ratio was also noted to be dif- 
ferent from that produced by either simple optical 
depth effects or chromospheric models (Hamann 
& Persson 1992). We have compared the obser- 
vations and model predictions presented in both 
Fig. 8 of Hamann & Persson (1992) and in Fig. 3 
of Herbig & Soderblom (1980) with our derived ra- 
tios to investigate the physical conditions in which 
Ca II lines arise. We note that in the above pa- 
pers, the quantities used are Wa (8498)/ Wa (8542) 
and Wa(8662)/Wa(8542). Our measured ratios 
of -1.7:2:1.9 therefore become 0.86±0.05 and 
0.97±0.05, respectively. A similar ratio was also 
found by Gibb et al. (2006) and Vacca et al. 
(2004) in their post-outburst spectra. In Fig. 8 
of Hamann & Persson (1992), V1647 Ori would 
lie close to the loci of models for either i) a very 
high optical depth (r '--^lOOO) isothermal LTE slab, 
or ii) a somewhat lower optical depth isothermal 
non-LTE slab (r '^50, but still optically thick). 
The location of V1647 Ori is distinctly separated 
from both the majority of T Tauri stars they ob- 
served (i.e. with mean ratios of —1,0.8), and chro- 
mospheric model values (i.e. 0.6,0.8). Following 
Hamann & Persson (1992), ratios close to unity as 
in our case, indicate that coUisional decay is dom- 
inant over radiative decay and requires conditions 
with rUe >10^'^ cm~'^ where ne is the electron 
density. It seems, therefore, that regardless of 
LTE or non-LTE conditions, in February 2007 the 
Ca II triplet lines in V1647 Ori were extremely 
saturated and produced in a very high density 
medium. This is also supported by the location 
of V1647 Ori in Fig. 3 of Herbig & Soderblom 
(1980), where V1647 Ori would he close to the 
optically thick LTE slab location. At this lo- 
cation, several T Tauri stars are found, specifi- 
cally, GW Ori (0.93,0.95), V380 Ori (0.88,0.87), 
V1121 Oph (0.85, 0.91), and RW Aur (1.02,0.84). 
We finally noting that it appears the Ca II triplet 
ratio changed sometime during the post-outburst 
phase between March and November 2004 suggest- 
ing a fundamental physical change in the region 
where the Ca II lines are emitted. 



3.3. The Near-IR Appearance 

Of primary interest in our H-band image shown 
in Fig. [3] is the fact that V1647 Ori itself is neb- 
ulous and not just a point-source. The nebu- 
losity takes the form of curving tails extending 
from the star to the north-east and north-west. 
These features were also observed in the K-band 
by Reipurth & Aspin (2004), in J, H, and K by 
Ojha et al. (2004) and in J and K by AP07. Our 
image shows somewhat more extensive nebulosity, 
particularly at very faint levels, reaching as far as 
HH 22-IRS. These curving 'cometary' tails most 
likely trace the cavity walls of a remnant outflow 
cavity (although, as we have noted, no such out- 
flow has yet been found). However, it has also 
been suggested that such a cavity may be a rem- 
nant of infall to the circumstellar envelope (Whit- 
ney & Hartmann 1993; Kenyon et al. 1993; and 
Hartmann, Calvet, & Boss 1996). 

Such cavities, be they from infall or outflow, 
have been modeled by Reipurth et al. (2000) and 
later Stark et al. (2006). The purely geometric 
models of Reipurth et al. (2000) were appHed 
to HST NICMOS images of Herbig-Haro energy 
sources and predicted the cavity shape for different 
opening angles and inclinations. A cursory com- 
parison between the McNeil's Nebula optical im- 
age from February 2004 (Fig. [T]) and their results 
suggests that an inclination with respect to the 
line-of-sight, 6*, of ~60 degrees and cavity open- 
ing angle, a, of ~30 degrees is appropriate. More 
quantitatively, applying formulae (1) for and (2) 
for a from Reipurth et al. (2000) we obtain values 
of ^65 degrees and a ~30 degrees. 

In the models of Stark et al. (2006), the cav- 
ity walls are assumed to be seen in reflected light. 
They predict the near-IR outflow cavity geome- 
try for a variety of physical parameters associ- 
ated with both the circumstellar disk and an ex- 
tended infalling envelope. They also consider the 
effect of cavity inclination and the shape/form of 
the cavity e.g. a 'streamline' cavity with straight 
walls extending away from the central star, or 
a 'curved cavity' (as we see in V1647 Ori) on 
the derived morphology. Although these models 
were created specifically to interpret HST NIC- 
MOS images of young stellar objects in Tau- 
rus (e.g. IRAS 04302+2247, DC Tau B, and 
CoKu Tau/1), their use as a qualitative compari- 
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son to our VI 647 Ori H-band image seems appro- 
priate. Of the NICMOS observations considered 
in Stark et al. (2006), those of CoKu Tau/1 ap- 
pear most similar in morphology to our near-IR 
data on V1647 Ori although the spatial scales con- 
sidered are clearly much different. The numerical 
model that best fits the CoKu Tau/1 data involves 
a curving (bipolar) cavity with opening angle 20 
degrees and inclination to the line-of-sight of 64 
degrees, an infalling envelope with infall rate of 
2x10"'^ M0yr~^, and a disk with a 50 AU radius 
disk. Clearly, we cannot make any claim that the 
best- fit parameters for CoKu Tau/1 are appropri- 
ate for V1647 Ori, only that such parameters pro- 
duce a model that qualitatively resembles the cav- 
ity walls observed. An important point, however, 
is that the models predict that the extent of the 
cavity wall illumination changes significantly with 
outflow axis inclination and that the most appro- 
priate model for CoKu Tau/1 has an inclination of 
64 degrees, a value close to that derived by AP07. 

In addition to V1647 Ori, the H-band im- 
age shown in Fig. [3] contains several other point- 
sources including the interesting near-IR object 
HH 22-IRS. The nebulous extensions seen in the 
IRAC images of this source are also detected in 
our H-band images and likely take the form of 
cavity walls extending to the west from HH 22- 
IRS. These cavity walls appear more streamlined 
than curved. Higher spatial resolution, high sig- 
nahnoise data may shed light on the origin of both 
types of cavities. 

3.4. The Near-IR Spectroscopic Features 

Fig.|3]shows the near-IR spectrum of V1647 Ori 
from 1 to 4.1 /im taken with SpeX on the NASA 
IRTF telescope. The top panel shows the com- 
plete spectrum while the middle and bottom pan- 
els show expanded views of the J and H band 
(middle), and K band (bottom). Numerous spec- 
tral emission and absorption features are identified 
from both atomic and molecular transitions. In- 
spection of Fig. m shows that a year after the out- 
burst had subsided, V1647 Ori shows broad wa- 
ter vapor absorption features (1.3-1.6 fj.m, 1.75- 
2.2 /xm, and 2.35-2.5 /itm), as well as deep wa- 
ter ice absorption (2.8-3.3 /Ltm). Particularly sig- 
nificant is the presence of weak CO overtone ab- 
sorption (2.294 /im), weak atomic Na I absorption 
(2.206 /im), and weak Ca I absorption (2.264 /im). 



These features are best seen in the expanded 2.1- 
2.4 /im spectrum in Fig. [5l In emission, we see 
relatively strong Br7 (2.166 /im), weak H2 v=l-0 
S(l) (2.122 /im), weak [Fe II] (1.275, 1.533, 1.644, 
and 1.809 /J.m), and several other H I lines from 
the Paschen and Brackett series, the longest in 
wavelength being Bra (4.05 /im). 

Of particular note is the evolution of the CO 
overtone features from 2.3-2.5 /im through the 
outburst. In the near-IR spectra of Reipurth & 
Aspin (2004) and Vacca et al. (2004), taken in 
February and March 2004, respectively, the CO 
overtone bandheads were strongly in emission. As 
the outburst progressed the CO emission band- 
heads faded and in the May 2006 spectrum of 
AP07 the 2.3-2.5 /im region was featureless and 
remained so in their August 2006 spectrum. In 
our spectrum from February 2007, however, we are 
now seeing weak CO bandheads in absorption. As 
we shall see below (Section 3.5), we believe that 
we are observing CO, Na I and Ca I absorption 
from the stellar photosphere of V1647 Ori. 

The Br7 emission flux from V1647 Ori was also 
tracked by AP07. They found an approximately 
linear decline in Br7 flux from 2.2x10^^® W m^^ 
in March 2004 to 3.3x10"" W m''^ in September 
2006. Our Br7 flux of 2.3x 10"" W m^^ in Febru- 
ary 2007 follows this trend. A similar trend is 
found in both He I (1.083 /im) and Pa^ (1.282 /im) 
emission with an approximately linear trend de- 
clining by factors ~ x30 and xl7, respectively, 
from March 2004 to February 2007. 

The broadest features present in the near-IR 
spectrum of V1647 Ori are due to water molecules. 
Specifically, we see water vapor absorption fea- 
tures in the J, H, and K bands, and water ice 
absorption in the L band. While the water ice ab- 
sorption feature has been observed in V1647 Ori 
from the very first observations at 3 /xm (Vacca et 
al. 2004), water vapor absorption is detected for 
the first time. 

We find that the water ice absorption depth 
and profile in our spectrum appears very similar 
to the profiles found in the spectra presented by 
Rettig et al. (2005) and Gibb et al. (2006). This 
perhaps suggests that the ice mantle grains were 
not significantly affected by the outburst itself. 

The presence of water vapor absorption sug- 
gests that we are observing a cool photosphere 
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of spectral type later-M. However, Shiba et al. 
(1993) suggested that a significant fraction of 
young low-mass stars show water vapor absorp- 
tion in excess of their true spectral type. This was 
interpreted as evidence for absorption in the ex- 
tended atmosphere of a circumstellar disk. Qual- 
itatively the water vapor absorption we observe 
in V1647 Ori seems inconsistent with the atomic 
lines and molecular band absorption seen in the K 
band. A more detailed discussion of this is given 
below where we apply stellar models to the ob- 
served spectrum to derive physical parameters as- 
sociated with the young star itself and the circum- 
stellar material. 

3.5. Near-IR spectral modeling 

The detection of neutral Na (2.206 fim) and Ca 
(2.265 fim) atomic absorption, and ro-vibrational 
molecular CO v=2-0 bandhead (2.294 /itm) absorp- 
tion in the near-IR spectrum of V1647 Ori can be 
interpreted as evidence that we are detecting flux 
from a stellar photosphere. However, it is also 
possible that these spectral features are created 
in a disk atmosphere. The atmospheres of cir- 
cumstellar disks are considered to be low-surface 
gravity environments similar in many ways to the 
atmospheres of cool supergiants. In such atmo- 
spheres, molecular CO bandhead absorption has 
been found to be extremely strong with respect 
to neutral atomic Na and Ca absorption. Exam- 
ples of such disk-dominated near-IR spectra can 
be found in FUors where the 2/im spectrum is 
often relatively featureless besides deep, broad- 
ened CO absorption (see for example Fig. 1 in 
Rcipurth & Aspin 1997). In V1647 Ori, unless 
the CO absorption is significantly diluted by CO 
emission, the strength of the Na and Ca lines with 
respect to the CO absorption suggests that the en- 
vironment of the absorption region is more dwarf- 
like than supergiant-like. However, we find it is 
not perfectly dwarf-like but intermediate between 
dwarf and giant values. This is consistent with 
V1647 Ori being a young, pre-main sequence dwarf 
star still contracting along Hayashi tracks. In the 
discussion and analysis that follows, therefore, we 
assume that the Na, Ca, and CO absorption fea- 
tures are from the stellar photosphere of V1647 Ori 
but note that this assumption is not unambigu- 
ously determined. 

By modeling the strength of photospheric ab- 



sorption features, using, for example, spectral 
template stars, allows us to derive estimates of the 
underlying stellar type, the characteristics of the 
IR excess emission, and the extinction toward the 
photosphere. We have first applied the modeling 
method of Prato, Greene & Simon (2003) to the K- 
band spectra of V1647 Ori. This is the spectral re- 
gion where the aforementioned photospheric lines 
are found. The spectrum was fit with a series of 
models of spectral type template stars which were 
veiled with infrared emission excess and reddened 
by extinction using the ISM extinction law: 

= A„[0.55/A]i-6 (1) 

taken from Prato, Greene & Simon (2003) 
where 0.55 and A are in microns. The K-band 
veiling, tk, is defined as the ratio of the infrared 
excess at 2.22 /im to the intrinsic stellar contimmm 
flux. The resuhs of MuzeroUe et al. (2003) and 
Beck (2007) show that IR excess can be described 
by blackbody emission with a single character- 
istic temperature and hence, for the purpose of 
the K-band models, we use an IR excess with a 
blackbody temperature, Tdust, of 1000 K. A 
search for the optimum tk value and the best-fit 
Ay was performed for a range of stellar template 
stars from K5 to M3. The reader is referred to 
Prato, Greene & Simon (2003) and Beck (2007) 
for more details of the models. 

In Fig. [71 we plot the K-band spectrum of 
V1647 Ori (in black) with the best-fit models for 
three template stellar types, K7 V, MO V and 
Ml V. Stellar types in the K5 to M2 range provide 
the best average fit to the detected photospheric 
features: stars with much earlier types (i.e. ~K3 
or earlier) have Mg I (2.281 fim) absorption which 
should have been detected at the same strength as 
Ca (2.265 /im), and stars with much later stellar 
types (i.e. ~M3 and later) have much stronger 
Na I (2.206 /xm) absorption relative to Ca I. This 
is true for both dwarf stars and giants. The stel- 
lar templates were taken from the NASA IRTF 
SpeX stellar spectral library located on the IRTF 
web siteH. From the K-band model fits, the best- 
fit values give A^, ^--^29 magnitudes, and rfc=1.4- 
1.5. However, early in the process of modeling 
the K-band spectra, it became apparent that the 



^http: / /irtfwcb. ifa.hawaii.edu/ spcx/spexlibrary/IRTFlibrary.html 
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spectrum at wavelengths greater than ~2.3 /zm de- 
viated significantly from the models which best-fit 
the photospheric absorption features (see Fig. [7]). 
This deviation is best explained by excess absorp- 
tion from water vapor often seen, for example, in 
FUors (Greene & Lada 1996). The presence of wa- 
ter vapor features significantly affects the extinc- 
tion values determined by modeling the K-band 
alone. Fitting the steep slope of the K-band spec- 
tra at A <2.3 /xm (made steeper by the 1.9-2.2 /xm 
water vapor band) requires extinction value of ~29 
magnitudes. This is a much larger Ky than pre- 
viously associated with V1647 Ori (i.e. 6-14 mag- 
nitudes, Aspin & Reipurth 2004; Briceho et al. 
2004; Abraham et al. 2004; Andrews et al. 2004; 
Vacca et al. 2004, Ojha et al. 2006, AP07) and, 
since water vapor absorption occurs shortward of 
2.2 ^m also, we conclude that the K-band models 
are likely affected by the presence of water vapor 
and are therefore inadequate for the derivation of 
physical characteristics. 

To further investigate the above models valid- 
ity, we have extended the near-IR spectral mod- 
eling of Prato, Greene, & Simon (2003) and the 
version used by Beck (2007) to include the full 
1.0-4.1 fiuY spectral range in the hope of placing 
more stringent constraints on both Ay and vk- 
The extended synthesis model was adapted from 
the K-band spectral code to use an optical depth 
model which included water ice absorption, and an 
infrared excess model that included water vapor 
absorption. The model of the observed 1-4.1 /im 
spectrum is defined as: 

Fobs = Cx [{F„u + fc(A))e-^(^)] (2) 

where Fi„t is the intrinsic photospheric flux of 
the young star estimated from the spectral tem- 
plate star. The wavelength dependent optical 
depth, t(A), has two components: first from the 
ISM extinction law as described in Eqn. 1 above, 
and second from an optical depth component, r^ce, 
derived directly from laboratory absorbance data 
for water ice (Gerakines et al. (1995; 199^1). The 
parameter C is a scaling constant to normalize the 
model to unity at 2.22 ^m. The value k(A) is the 
wavelength dependent emission excess, which is 
described using black-body radiation with a tem- 



Laboratory water ice absorbance data are available online 
at: |http: / /www.strw.leidenuniv.nl/^lab| 



perature, T^ust, constrained by the required tk 
value at 2.22 /im determined by the K-band mod- 
els (Muzerolle et al. 2003). We have made an 
important modification to the Muzerolle et al. in- 
frared excess model by incorporated an estimate 
of water vapor absorption. To accomplish this, 
we have multiplied the blackbody curve that de- 
scribes the infrared excess emission by a smoothed 
template model of water vapor absorption. The re- 
sults are presented in Fig. [S] The red curve is the 
blackbody IR excess with a temperature of 1000 K, 
and the blue curve is this blackbody multiplied 
by the water vapor absorption model. The water 
vapor absorption model was constructed hy i) re- 
moving all photospheric absorption features from 
a M7 V star (using the SpexTool "xcleanspec" 
routine - Gushing, Vacca, & Rayner 2004), ii) 
smoothing the cleaned spectrum using a 10-pixel 
boxcar. Hi) dividing the smoothed spectrum by a 
black-body of temperature 2600 K to remove the 
M7 V spectral continuum shape, and iv) normal- 
izing the flux of the water vapor model to unity at 
2.22 /im. An M7 V stellar spectrum was adopted 
since it best fit the water absorption observed in 
V1647 Ori. Since water vapor characteristics are 
not derived by the fitting procedure and only al- 
low us to obtain more meaningful values of the 
free-parameters, we deemed this methodology ac- 
ceptable. 

A minimization search was performed over 
the three parameters, A^, Tdust and r^ce as in Beck 
(2007) . To obtain meaningful uncertainties on the 
derived values, we have generated a range of 
surfaces for A„ and Tdust values for the K5 to 
M2 range of stellar types which best fit the K- 
band spectrum (e.g. Fig. [7]). A minimization 
search was performed over the three parameters, 
A«, Tdust and Tice- As found by Beck (2007), the 
optical depth of the water-ice absorption feature 
did not affect the models appreciably in that only 
small changes in r^ce are seen for a wide range of 
A^, and Tdust values. For the analysis involving 
a range of stellar types, the infrared veiling value 
was adopted to be rfc=1.5 and the shape of the 
IR excess incorporated water vapor excess as in 
Fig. [H Fig. [9] presents the x^ = l-01, 2.0 and 3.0 
surfaces in A^ and Tdust for the models using K5 V 
(red), K7 V (green), MO V (blue). Ml V (cyan), 
and M2 V (magenta) spectral template stars. Be- 
cause the models are fitting a broad spectral range. 
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the Icr surface (approximated by x^=3.0) is very 
narrow in Ay and Tdust- Best-fit At, values range 
from 17.5 to 21.5 magnitudes, while best-fit Tdust 
values span the range 900 K to 1150 K. Based 
on these results and the K-band fits, we adopt the 
MO V stellar template star to be the most accurate 
model, with an uncertainty of ±2 spectral classes. 

To further define the uncertainties in the spec- 
tral models, a second set of searches was per- 
formed over A„ and Tdust- For this search, a spec- 
tral type of MO V was adopted and the infrared 
veiling value, r^, was the free parameter. Fig. 1101 
presents the x^=1.01, 2.0 and 3.0 surfaces for 
vs. Tdust using rfe=1.3 (red), 1.4 (green), 1.5 
(blue), 1.6 (cyan) and 1.7 (magenta). Again, for a 
range of infrared veiling values, models that fit the 
broad 1.0-4.1 /im spectral range show very small 
variation in the surface. Based on these analy- 
ses, the best- fit model for the 1.0-4.1 /im near-IR 
spectra of V1647 Ori (including the effects of wa- 
ter vapor absorption) is adopted to be an M0±2 
star with Au=19±l magnitude, rfc=1.5±0.2, and 
a Td„^f=1000 KiblOO K. 

Fig. [TT] presents the full 1.0-4.1 fj,m spectrum 
of V1647 Ori with the best-fit model overplot- 
ted in red. Unfortunately, the models cannot ex- 
plain the known long-wavelength absorption wing 
to the 3.0 /im water-ice absorption feature (the 
3.27-3.70 fim region was excluded from the spec- 
tral model fits - Smith et al. 1989; Beck 2007). 
Fig. [H shows the fuU x^ = l-01, 2.0, 3.0, 6.0, and 
10.0 surfaces for the model using an MO V tem- 
plate star. The x^=10.0 surface is a measure of 
the 3(7 uncertainties in the model. Visual extinc- 
tions as low as '^11 magnitudes A^, as measured by 
other studies (e.g. AP07) are excluded at a >3ct 
level of confidence in our models. 

3.6. Additional determinations of visual 
extinction 

Of crucial importance to the determination of 
the physical parameters of V1647 Ori is the deter- 
mination of reliable visual extinction values. With 
the dramatic optical fading of V1647 Ori of over 
5 magnitudes. Ay values towards the young star 
likely vary with time as suggested by Rcipurth & 
Aspin (2004). McGehee et al. (2004) determined 
that at optical wavelengths the change in color of 
the source was a result of both an Ay change and 
the addition of significant blue accretion luminos- 



ity. However, in the near-IR, the effect of such 
an increase in optical brightness would probably 
be much smaller. Below, we examine our near-IR 
photometry and spectroscopy and attempt to de- 
rive Ay values using several different and indepen- 
dent methods and compare them to Ay ^19±1 
derived from the spectral template fitting. 

3.6.1. Ay from Near-IR colors 

From our calibrated near-IR spectrum we can 
obtain estimates of the broad-band J, H, and K 
magnitudes of V1647 Ori. The determination of 
J and K magnitudes from this spectrum was per- 
formed using our ULBCam H band imaging mag- 
nitude of 11.96 to determine an absolute calibra- 
tion. Using the conversion defined by Tokunaga & 
Vacca (2005) for J and K bands, we obtain J=14.8, 
H=11.96, and K=9.87 and hence J-H=2.8 and H- 
K=2.1. Dereddening these values along a redden- 
ing vector to the CTTS locus defined by Meyer et 
al. (1997) gives an Ay ~17±4 magnitudes where 
the uncertainty is derived from the errors on the 
photometric values. 

An interesting feature of the published near- 
IR colors for V1647 Ori is that values from Octo- 
ber 1998 (2MASS), March 2004 (Reipurth & As- 
pin 2004), and May 2006 (AP07) all he close to 
the same reddening vector and have a range of 
Ay values from ~9 (at the peak of the outburst 
phase) to ~17 magnitudes (one year into the qui- 
escent phase). This trend demonstrates that the 
Ay towards the star decreases during outburst 
and increase as the source fades toward quies- 
cence. Fig. ll3l shows this graphically where we plot 
published J-H and H-K colors covering the pre- 
outburst, outburst and the quiescent phase (our 
data). Using these published near-IR colors, we 
have derived the Ay values for each observation by 
dereddening onto the CTTS locus. The estimated 
Ay values are plotted against time in Fig. [TH 
It is clear from this plot that the Ay towards 
V1647 Ori is considerably larger in its quiescent 
phase than during the outburst phase. We have 
considered whether a change in Tdust near the star 
could mimic the change in Ay seen in Fig. 1141 
Near-IR J-H and H-K colors for pure Planck func- 
tions with temperatures in the range 950 K to 
1350 K do give similar colors as V1647 Ori exhibits 
in Fig. [ill however, for this to be the underly- 
ing cause of the observed color changes from out- 
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burst to quiescence wc would have to be observing 
purely thermal emission from the cool black-body 
with an insignificant contribution from the stellar 
photosphere. Since we detected photospheric fea- 
tures in the quiescent phase K-band spectrum of 
V1647 Ori, this interpretation is clearly invalid. 

It also appears that i) the near-IR colors of 
V1647 Ori are little effected by the addition of the 
hot accretion component to the luminosity during 
the outburst phase, and ii) that the near-IR ther- 
mal excess remained approximately constant from 
the outburst to the quiescent phase with a change 
in visual extinction of AAy ^8 magnitudes. The 
only way (i) above can occur is if the accretion lu- 
minosity is 'gray' at near-IR wavelengths i.e. the 
accretion luminosity increases the emergent flux 
but does not change the observed near-IR colors 
of the source. Gray absorption has been previously 
observed in accreting young stars in Orion A by 
Carpenter, Hillenbrand, & Skrutskie (2001) and 
interpreted as being due to large particles. How- 
ever, whether a similar mechanism can be invoked 
for emission is unclear. As to (ii) above, it would 
seem that the level of accretion in the pre-outburst 
phase (when the near-IR colors were measured by 
2MASS i.e. October 1998) must have been similar 
to that in quiescence (in February 2007), otherwise 
the SEDs would not be close to identical. 

3.6.2. Ay from near-IR [Fe II] line ratios 

Reipurth et al. (2000) used the two near-IR 
[Fe II] lines at 1.275 and 1.644 /im to obtain an es- 
timate of the visual extinction towards the [Fe II] 
emitting region in a number of HH energy sources. 
This is possible since the transitions forming the 
lines arise from the same upper level and therefore 
the line ratio is fixed. We note that the estimate 
obtained will be the line-of-sight extinction to the 
[Fe II] emitting region which may be located some 
distance from the stellar photosphere and possi- 
bly extended in nature. Nonetheless, using the 
formula shown in Reipurth et al. (2000), inverted 
to estimate Ay i.e. 

Av = -25 X Zoffio((/l.275//l.644)/1.36) (3) 

for an extinction law that gives Ay ~10xEj_h, 
the [Fe II] line fluxes in Table 2 give a ratio 
Ii. 275 /Ii. 644 =0.3 and therefore an Ay '~16±3 mag- 



nitudes. The error is estimated from repeat mea- 
surement of the (weak) line fluxes. 

3.6.3. Ay from the 3 [im water ice absorption 

feature 

As outlined in Beck (2007), deriving Ay val- 
ues from the 3 /zm water ice absorption feature 
is somewhat unreliable for Ay >20 magnitudes 
due to uncertainties in the relationship between ice 
band optical depth and Ay. Also, grains nearer 
to the star that still contribute to the line-of-sight 
extinction, may not have ice mantles if the local 
temperature is >200 K. Nonetheless, since all in- 
dications arc that for V1647 Ori we are dealing 
with Ay values <20 magnitudes (e.g. AP07), we 
proceed with the derivation. 

Following Beck (2007) and using the relation 
presented in Teixeira & Emerson (1999) between 
water column density, N(H20), and Ay, for a 
peak optical depth in the water band Tp^ak ~0.7, 
and a FWHM of the band of 320 cm.-\ we ob- 
tain N(H2O)~1.0 10^* cm-2. This value is almost 
identical to that derived by Rettig et al. (2006) i.e. 
1.2 10^® cm~^. From Fig. 5a in Teixeira & Emer- 
son (1999), we see that the above column density 
corresponds to an Ay '^13±4 magnitudes. How- 
ever, as found in Beck (2007), changing the value 
of Ay used resulted in only small changes in the 
model continuum level and did not substantially 
change the derived ice absorption optical depth. 
In light of this, we consider that the Ay derived 
using the ice absorption optical depth is likely not 
as accurate as the values from the other determi- 
nations above. 

3.6.4- from near-IR H I line ratios 

Another method for determining Ay towards 
the emitting region is the measured ratio of the 
H I lines Br7 (2.166 /im) and Bra (4.052 /im). 
This determination requires the assumption that 
only Lya photons are optically thick (i.e. Case B 
of Baker & Menzel 1938) over a wide range of tem- 
peratures and densities. If that assumption is cor- 
rect then the ratio Bra/Br7 should be in the range 
2.85-3.17 (Hummer & Storey 1987). Due to the 
relatively low signal to noise of Bra it is diSicult to 
estimate the line flux accurately. Our best value of 
the above ratio, obtained from repeated measure- 
ments of the Bra line flux, is 2.0±0.6. This sug- 



13 



gests that the Bra hne flux is less that expected 
from Case B and, hence, the Case B assumption 
is Ukely invahd. As a result of this, we do not at- 
tempt to derive Ay values from the H I emission 
ratios. 

3.6.5. Summary of derived Ay values 

To summarize, our four independent determi- 
nations of Ay i.e. from spectral template fit- 
ting (19±1), near-IR colors (17±4), [Fe 11] hne 
ratios (16±3), and 3 /xm ice absorption optical 
depth (13±4), all give reasonably consistent re- 
sults. However, since we put lower weight on 
the value determined from the 3 /xm ice absorp- 
tion feature, and are most confident in the Ay 
value derived from our spectral template fitting, 
we adopt Ay ~19±2 magnitudes as the visual ex- 
tinction towards V1647 Ori in February 2007 and 
note that the Ay values derived from near-IR col- 
ors and [Fe II] line ratios are consistent with this 
value (within the associated uncertainties). Addi- 
tionally we have found that there was a systematic 
decrease in Ay towards V1647 Ori during the out- 
burst phase. This was mirrored by an increase in 
Ay as the source faded to, and into, quiescence. 

3.7. Derived physical parameters 

3.7.1. Luminosity and absolute magnitude 

Above, we determined that V1647 Ori is likely 
a star of spectral type M0±2 with a Teff of 
~3800 K, a visual extinction Ay ~19±2 magni- 
tudes (in February 2007), a veiling of rx=1.5±0.2, 
and an apparent K magnitude of m^ ^9.90 (in 
February 2007). Using the bolometric correction 
(M6oi-M„=-1.3) and MO color (V-K=3.7) from 
Table A5 in Kenyon & Hartmann (1995) and the 
equations for Mboi and from Greene & Lada 
(1997), we estimate that V1647 Ori has an abso- 
lute bolometric magnitude of M^oi '^2.9±0.4 and 
a luminosity of ~5.2±2 Lq for a distance of 
d=450 pc (m-M=8.3). The error on our value of 
L* takes into account the uncertainties on tk , Ay , 
and the spectral type quoted above, but not un- 
certainties in distance. Thus, we estimate that the 
star has a radius, R, , of ^5 Rq using the standard 
Te// vs. relationship. 



3.7.2. Stellar mass and age 

Using the stellar luminosity and spectral type 
derived above, we can place V1647 Ori on pre- 
main sequence stellar tracks using results from 
numerical models. Our adopted values are 
L* ~5.2±2 Lq and Te// -3800±350 K. The 
temperature (and associated uncertainty) comes 
from our adopted spectral type of MO (±2 sub- 
classes) using the values in Table A5 of Kenyon 
& Hartmann (1995). We have placed V1647 Ori 
on the evolutionary tracks of both Baraffe et al. 
(1998) and Siess et al. (2000). Figs. HH and [H 
show its location on these model tracks. The 
models give masses and ages for V1647 Ori of 

1 Mo/<0.5 Myrs, and 0.55 Mq/<0.5 Myrs, re- 
spectively. We conclude therefore that, conserva- 
tively, V1647 Ori is probably of mass 0.8±0.2 Mq 
with an age of < 0.5 Myrs. 

3.7.3. Mass accretion rate and accretion lumi- 
nosity 

It is possible to determine mass accretion rates 
for young stars like V1647 Ori using H I emis- 
sion line strengths. Following MuzeroUe, Hart- 
mann, & Calvet (1998), we can use the dereddened 
fiux of the near-IR emission lines Pa/3 and Br7 to 
derive accretion luminosities, haca and following 
GuUbring et al. (1998) determine a mass accre- 
tion rate, M. The reader is referred to Eqns. 1 & 

2 in MuzeroUe, Hartmann, & Calvet (1998) and 
Eqn. 8 in GuUbring et al. (1998) for further de- 
tails of the derivation. For the calculations below, 
we assume Mstar=0.8±0.2 Mq, L*=5.2±2 L©, 
T,t„^=3800 K, R„,=5 Rstar, Ay=19, the extinc- 
tion law shown in Eqn. 1 above, and d=450 pc. 
Ri„ is the inner radius of the accretion disk. 

Dereddened fiuxes and derived Ax values for 
all detected optical and near-IR emission lines are 
shown in Table 2. Using the dereddened Pa/3 
fiux, we obtain L^cc '~~^3.5 Lq and log(M) — 6.1 
or M~8.8 lO"'^ Mq yr'^ With the uncertain- 
ties associated with the parameters in the for- 
mula used, the range of possible values for L^cc 
are 2.1-5.5 Lq and for M is 5.4 lO"'' to 1.4 lO"*^ 
Mq yr~^. Similarly, for the dereddened Br7 fiux, 
we obtain Lace ~4.6 Lq and a log(M)~-5.92 or 
M~1.2 10-6 Mq yr-i. Considering the quoted 
uncertainties, we obtain a range of possible val- 
ues of Lace of 2.8-7.7 Lq and M=7.2 10"'^ to 
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2.0 10-6 M© yr-i. 

Another method to derive M was discussed by 
White & Basri (2003) and Natta et aL (2004). 
This derivation relates the width of the optical 
Ha emission line to accretion rate. The quan- 
tity used for this is the Ha full width at 10% 
peak intensity and an approximately linear rela- 
tionship was found between this and log(M) by 
Natta et al. (2004). Even though the Ha hue pro- 
file appears asymmetric we apply this technique 
to investigate its applicability. Using Eqn. 1 from 
White & Basri (2003), and a measured Ha at 10% 
width of 740 km s"\ we obtain a log(M)~-5.7 
or M~2 IQ-^ M0 yr"^. This derivation assumes 
that the Ha emission is solely the result of the ac- 
cretion processes and not contaminated by stellar 
winds (absorption/emission) and/or HH flows/jets 
(emission) which, due to the presence of an asym- 
metric profile, may not be the case. However, we 
note that if we double the blue-wing 10% emission 
width to compensate for the red-shifted absorp- 
tion, we would obtain a larger value of Ha at 10% 
and hence, a larger accretion rate. 

To summarize, our two estimates of accretion 
luminosity and rate give reasonably consistent re- 
sults and we cldopt VcllueS of Lace 

--4±2 L0 and 
M~1.0±0.5 10-6 M© yr-1, respectively. 

3.7.4- Bolometric Luminosity 

Our estimate of hhoi is the sum of our derived 
L, and hacc values, or, L;,oi=9.2±3 L©. An- 
drew, Rothberg, & Simon (2004) calculated an 
outburst luminosity of Lf,oi '^3.5 L© and noted 
it was similar to the one derived by Lis, Menten, 
& Zylka (1999) once a distance correction was ap- 
plied (they assumed d=400 pc). They also con- 
sidered their luminosity estimate a lower limit due 
to the lack of information on the spectral energy 
distribution (SED) peak in the far-IR and the lin- 
ear interpolation used to define the SED. We note 
that the Lboj estimate from Lis, Menten, & Zylka 
(1999) only considered emission from two modified 
Planck functions that fitted the IRAS and sub- 
mm fiuxes and did not include a stellar compo- 
nent. In addition, it is clear that the IRAS fluxes 
can contain contributions from nebular emission 
and from other stars in the beam thus contami- 
nating the measurements. We also note that An- 
drews, Rothberg, & Simon (2004) did not quote 
any formal uncertainty on their hboi determina- 



tion. The error value is probably significant in 
light of the points related above. We consider our 
estimate of Ijboi to be robust in that it has consid- 
ered both stellar luminosity and accretion lumi- 
nosity independent of the source SED. Since the 
SED of V1647 Ori is very similar pre-outburst and 
in the quiescent phase (see Section 3.7.5), we con- 
sider that the Lt,oi of the source during both these 
periods is 9.2±3 L©. 

3. 7. 5. Spectral Energy Distribution 

The SED of V1647 Ori is shown in Fig. [ITl 
The open stars are the data from February 2007 
while the open boxes and solid points are the pre- 
outburst and outburst SEDs, respectively. We 
have used the same data as Andrews, Rothberg, & 
Simon (2004) longward of 25 jim (i.e. IRAS, and 
JCMT data) and note that they state that the 
pre-outburst and outburst photometry from 350- 
1300 fim was unchanged by the eruption. It is ob- 
vious that the quiescent phase SED is almost iden- 
tical in shape and flux level to the pre-outburst 
SED. This suggests that the state of V1647 Ori 
in February 2007 was very similar to that in the 
pre-outburst period. 

3.8. The Mid-IR Characteristics 

In Fig.[TH]we show the mid-IR N-band spectrum 
of V1647 Ori from UT 5 April 2007. These data 
were calibrated using the mid-IR photometric data 
presented in Table 2. The spectrum is featureless 
apart from a change in slope around ^^9.5 /im. 
Shortward of this wavelength the spectrum is rel- 
atively flat. At ^^9.5 /im, however, the slope in- 
creases rapidly to the edge of the band. A previous 
mid-IR spectrum from Andrews, Rothberg, & Si- 
mon (2004) taken in March 2004 did not show this 
change in spectral index. Their spectrum rose to 
the red with a power-law fit F,y ~ A^'^. A similar 
fit to the 9.2-13 fim region of our spectrum gives 
F, ~A2-8. 

One difference between our N-band spectrum 
and that of Andrews et al. is that the flux level has 
decreased significantly from ~7 Jy to ~0.14 Jy at 
10 /um, a decline of a factor x50. In a study of the 
pre-outburst infrared characteristics of V1647 Ori, 
Abraham et al. (2004) presented ISOCAM pho- 
tometry at 6.7 and 14.3 fiui. If we interpolate 
these values to 10 fim we obtain a pre-outburst 
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brightness of ~0.4 Jy. In addition, the IRAS 
12 iim point-source flux for IRAS 05426-0007 is 
0.53 Jy. Given the uncertainties in the nature of 
the emission included in both the ISOCAM and 
IRAS beams, we consider this a reasonable agree- 
ment. 

We consider that the aforementioned change in 
spectral slope at ~9.5 /im is in fact indicative of 
the presence of weak silicate absorption through 
the passband. We follow the analysis of Quanz 
et al. (2007), who studied a number of FUors 
in the mid-IR using ISO and IRS spectroscopic 
data. They fitted a linear continuum between 
8 /im and 13 /^m and subtracted the observed 
spectrum from the result of the fitting procedure. 
In Fig. [18] we show such a linear fit as a dashed 
line. Subtracting a dereddened (with Av=19) ob- 
served spectrum from this fit and converting the 
result from fiux to optical depth, r, we obtain the 
plot shown in Fig. [12] where a relatively shallow 
yet typical silicate absorption band seen with a 
maximum optical depth of t ^0.35. Converting 
this value to an extinction at the absorption min- 
imum i.e. 9.2 /im (wc termed this A9.2) results 
in a value of 0.38±0.15 magnitudes. For com- 
parison, the expected extinction at 9.2 /im for an 
Ay=19 is A9.2 ~0.2. Additionally, comparing our 
silicate absorption plot with those from Quanz et 
al. (2007, their Fig. 15) we see a number of FUors 
with similarly shaped band, specifically, V346 Nor 
and Z CMa, which also have similar maximum op- 
tical depth to V1647 Ori (of r --0.3). A cursory 
comparison of the silicate absorption band shape 
with the model silicate bands from Quanz et al. 
(2007, their Fig. 16) suggests that the V1647 Ori 
silicate band is most similar to that modeled using 
astronomical silicates as defined by Weingartner & 
Draine (2001) and Draine (2003). 

In the observations of Andrews, Rothberg, & 
Simon (2004) from UT March 11 2004, there is no 
evidence for any silicate feature be it in absorption 
or emission. The mid-IR spectra of Quanz et al. 
(2007), from three epochs (UT October 20 2004, 
UT March 11 2005, and UT March 24 2005) were 
all found to show weak 10 /im silicate emission 
features. Mosoni et al. (2005) shows an mid-IR 
spectrum of V1647 Ori taken UT March 1 2005. 
This showed a similar weak silicate emission fea- 
ture seen by Quanz et al. Fedele et al. (2007a) also 
presented 10 /im spectra from three epochs (UT 



March 8 2004, UT February 21 2005, and UT Jan- 
uary 10 2006). In these, we believe that silicate 
emission is weakly present in only the February 
2005 spectrum. The progression therefore seems 
to be that i) silicate emission developed from a 
featureless red continuum between UT March 11 
2004 and UT October 20 2004, ii) the emission 
persisted until at least UT March 24 2005, and in) 
by UT January 10 2006 there was perhaps weak 
silicate absorption which grew a little deeper as 
time progressed to the level shown in our 10 /^m 
spectrum from UT April 05 2007. This result im- 
plies that there must have been significant dust 
evolution in V1647 Ori over the outburst period 
from March 2004 to April 2007. 

It is well known that some young eruptive vari- 
able stars show a 10 /im silicate feature in emission 
e.g. the FUors FU Ori, BBW 76, and V1057 Cyg 
(Green et al. 2006). Others, however, show 
a 10 /im silicate feature in absorption e.g. the 
FUors V346 Nor, Z CMa, L1551 IRS5 (Green 
et al. 2006). Silicate absorption is clearly pro- 
duced when dense, cool silicate dust is located 
between the radiation source and the observer. 
Silicate emission occurs when heated, optically 
thin, silicate dust, at a higher temperature than 
the source of the 10 /xm continuum emission, is 
directly visible to the observer. In addition to 
V1647 Ori, Quanz et al. (2007), studied of a num- 
ber of FUors, using Spitzer and ISO data, and 
suggested that two categories exist, one with ab- 
sorption created in a dense, dusty, and icy enve- 
lope (Category 1, herein CI), and the other with 
emission from the heated surface layers of a cir- 
cumstellar/accretion disk (Category 2, herein C2). 
They considered the progression from CI to C2 
evolutionary in nature and related to the size of 
the cold, icy dust envelope. Hence, CI FUors 
would be younger objects that still possess ex- 
tensive, dense, cold envelopes through which we 
are observing the FUor. C2 FUors would be more 
evolved in nature with only remnant envelopes and 
a more exposed accretion disk. Since FUor erup- 
tions are thought to repeat many times through 
the Class I to Class II stages, this is clearly a possi- 
bility although system orientation with respect to 
the line-of-sight may be important. Even though 
V1647 Ori is not considered a FUor, perhaps the 
clearing/sublimation of dust towards the star re- 
sulted in the accretion disk being revealed (C2) 
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and as the eruption and wind subsided significant 
dust re-accumulated/condensed along the hne-of- 
sight, obscuring the accretion disk, and creating 
enough column density to create the weak silicate 
absorption feature (CI). 

We defer further consideration of the 10 /um 
spectra to a future paper in which we will con- 
sider dust evolution more comprehensively using 
all published and archival 10 /xm spectra together 
with new 10 fim spectra taken during our two 
year Gemini monitoring campaign (Aspin et al., 
in prep). We conclude by noting thatzj the ap- 
pearance (and disappearance) of silicate emission 
supports the interpretation that the mid-IR flux 
originates from the circumstellar disk as suggested 
by MuzeroUe et al. (2005) and Abraham et al. 
(2006), and ii) the appearance of silicate absorp- 
tion suggests that the line-of-sight column density 
of cold dust increased as the outburst faded and 
the source returned to its quiescent state. 

4. CONCLUSIONS 

Prom our optical, near-IR, and mid-IR imaging 

and spectroscopy of V1647 Ori dating from Febru- 
ary to April 2007, approximately one year after its 
return to its pre-out burst optical brightness, we 
can conclude that: 

• The associated nebula, McNeil's Nebula, re- 
mained faintly visible suggesting that quasi- 
static nebula material is still being illumi- 
nated by V1647 Ori. 

• We confirm the findings of Fedele et al. 

(2007a) in that signposts of shock-excited 
emission, specifically, in the emission lines 
of [S II], [Fe II], and H2 are present. We 
consider that it is a distinct possibility 
that these emission lines result from a new 
Herbig-Haro flow. 

• Our near-IR spectrum shows, for the first 

time, molecular overtone absorption from 
CO and atomic absorption from neutral Na 
and Ca. We interpret this as evidence that 
we arc now observing the stellar photosphere 
of V1647 Ori. 

• We have modeled the near-IR spectrum us- 
ing a template stellar spectrum including 



both water vapor and water ice absorp- 
tion and have determined a best-fit param- 
eter set of spectral type M0±2 sub-classes, 
Ay=19±2 magnitudes, rif=1.5±0.2. 

• We derive values of Mf,o;=2.9±0.4 magni- 
tudes and L*=5.2±2 Lq for V1647 Ori. 
From comparison with theoretical evolu- 
tionary tracks, find that, for the adopted 
Teff ~3800 K, the star has a mass and age 
of 0.8±0.2 M0 and <0.5 Myrs, respectively. 

• From near-IR H I line flux and the Ay 
values derived above, we estimate the ac- 
cretion luminosity and mass accretion rate 
in February 2007 was Lacc=4.0±2 Lq and 
M= 1.0x10"*^ Mq yr-^ respectively. This 
implies that V1647 Ori is still actively ac- 
creting circumstellar material even though 
it is almost optically invisible and that the 
accretion rate during the outburst must have 
been considerably larger that this value. 

• V1647 Ori is found to have a quiescent phase 
Uoi of 9.25±3 L0. 

• For the first time, we see evidence for sili- 
cate dust evolution in the mid-IR spectrum 
of V1647 Ori over the outburst to quiescence 
period. We now observe weak silicate ab- 
sorption at 10 /xm whereas previously the 
silicate band was either absent or weakly in 
emission. 

• Finally wc note that, in February 2007, 
the spectral energy distribution, SED, of 
VI 647 Ori appears remarkably similar to its 
pre-outburst SED suggesting that perhaps 
the derived accretion rate is the normal qui- 
escent phase accretion rate for this object. 
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Table 1 
VI 647 Ori Observing Log 



UT Date 


MJD"* 


Telescope/Inst. 


Filters/Grism/Ac/Slit Width 


Exposure Times 


Seeing 


Photometry 


20051010/11 


3654/5 


Gcmini-S/GMOS-S 


[S II] 


1200 s 


0.65" 


- 


20070213 


4144 


UH2.2/ULBCam 


H 


360 s 


0.80" 


11.96'"±0.1'" 


20070221 


4152 


Gemini-N/GMOS-N 


R400/700 nm/0.75" 


3.5 hrs 


0.76" 


- 


20070222 


4153 


Gcmini-N/GMOS-N 


g' 


600 s 


0.66" 


26.21"±0.2™ 


zOOTOzzz 


4153 


Gemini-iN/GMUb-iN 


r 


OUU S 


U.53 


z3. ZD ±0.15 


20070222 


4153 


Gemini-N/GMOS-N 


i' 


600 S 


0.59" 


21.14™±0.1™ 


20070222 


4153 


Gemini-N/GMOS-N 


z' 


600 s 


0.63" 


19.05™±0.1™ 


20070222 


4153 


IRTF/SpeX 


SXD,LXD/0.8" 


1600 s,1200 s 


0.90" 




20070318 


4177 


Gemini-N /Michelle 


N'(11.2 Aim) 


100 s 


0.48" 


0.23±0.02 Jy 


20070318 


4177 


Gemini-N/Michcllc 


Qa(18.5 /im) 


100 s 


0.66" 


0.44±0.05 Jy 


20070405 




Gemini-N/Michelle 


lowN/0.4" 


800 s 


0.80" 





^Modified JuHan Date. 2450000-F 



Table 2 

Optical and Near-IR Equivalent Widths & Line Fluxes 



Line/Band 


Wavelength 


(A) 


— 

Observed Line Flux 


Dereddened Line Flux*^ 




( Mm) 


(W 


(W 


O I 


0.5579 


-104.4 


4.7 10 


1.3 10 


Ha 


0.6563 


-128.4 


1 10^^** 


5.4 10~^^ 


[S II] 


0.6716 


-16.3 


1.8 10-1" 


6.0 10"" 


[S II] 


0.6731 


-22.0 


2.3 10"^^ 


7.3 10"" 


Fe I 


0.7912 


-7.2 


4.0 10^^^ 


7.1 10"^'' 


O I 


0.8449 


-5.0 


3.8 10~^^ 


2.5 10"^'' 


Fc I 


0.8467 


-1.9 


1.7 10"^^ 


1.1 IQ-^^ 


Ca II 


0.8498 


-33.3 


2.9 lO""^* 


1.8 10"" 


Ca II 


0.8542 


-37.5 


3.4 10^^* 


1.9 10"" 


Fe II 


0.8600 


-2.7 


2.6 10""^^ 


1.4 10"^^ 


[Fc II] 


0.8617 


-3.0 


2.7 10"^'' 


1.4 10"^^ 


Ca II 


0.8662 


-35.1 


3.4 10-^® 


1.6 10~" 


Fe I 


0.8689 


-1.6 


1.6 10~^^ 


7.3 10"^'' 


Pa 12 


0.8748 


-3.1 


3.2 10~^^ 


1.3 10"^^ 


Pa 11 


0.8860 


-4.2 


4.6 10""-^ 


1.6 10"^^ 


Pa 10 


0.9012 


-4.7 


4.9 10""^^ 


1.4 10"^^ 


He I 


1.0830 


-10.0 


2.0 10~^* 


7.6 lO"^** 


Pa7 


1.0941 


-23.0 


3.5 10-^* 


1.2 10-^^ 


[Fe II] 


1.2570 


-2.0 


7.0 10~^^ 


7.4 10~" 


Pa/3 


1.2822 


-25.0 


8.6 10~^® 


7.9 10~^° 


[Fe II] 


1.5335 


-2.0 


2.3 10-^® 


6.8 10-1'' 


[Fe II] 


1.6445 


-2.0 


2.8 10-^®'^ 


5.8 10-" 


Br 10-4 


1.7367 


-3.0 


6.2 10-^® 


1.0 10-^** 


BiS 


1.9451 


-3.0 


1.0 10"^'' 


1.0 10-^** 


v=l-0 S(l) 


2.1218 


-1.0 


5.2 10-^* 


3.9 10-" 


Br7 


2.1661 


-5.0 


2.3 10-^'' 


1.6 10-^" 


Na I 


2.2060 


0.58 


absorption 


absorption 


Ca I 


2.2650 


0.49 


absorption 


absorption 


v=2-0 CO 


2.2935 


1.90 


absorption 


absorption 


Bra 


4.0522 


-12.0 


7.1 10-^'' 


1.5 10-^" 



^Wa is the line equivalent width. Negative Wa are emission, positive W>. are absorption. 

'^Emission line fiuxcs. Uncertainties arc estimated at ~ ±10% 

'^Emission line fluxes dereddened by Ay = ^9^ using AA=Av(0.55/A(//m))^'^ 

'^Possible inverse P Cygni profile. 
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Fig. 1. — Two images of V1647 Ori and McNeil's Nebula. The left image is a composite of g', r', and i' 
images taken on UT February 14, 2004. The right image is a composite of r', i', and z' images was taken on 
UT February 22 2007. North is at the top, East to the left in both images. The images are approximately 
70" by 80" in size at a pixel scale of 0.144"/pixel. For comparison, V1647 Ori, at the apex of the nebula, is 
~5.6 magnitudes brighter in the Feb 2004 image taken a few months after the initial outburst. 
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Fig. 2.— Optical spectroscopy of V1647 Ori taken on Gemini North using GMOS UT 22 February 2007. 
The R400 grating was used with a Ac=7000 A. This resulted in a spectral resolution, R, of ~1280. The top 
panel shows the complete optical spectrum of V1647 Ori. Note, the turnover in flux at 8900 A is an artifact 
of the data reduction. The middle panel is a zoom of the region around Ha, while the bottom panel is a 
zoom around the Ca II triplet lines. Flux, F^, is in W m~^ /im"^ and has been normalized to unity at 
8000 A by division by 5.1 xlQ-^^. 
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R.A. (J2000) 



Fig. 3. — A near-IR H band image of the region containing V1647 Ori and McNeil's Nebula. This image 
was taken on the University of Hawaii 2.2m telescope on UT 13 February 2007. V1647 Ori, the HH 22-IRS 
source, and Lk Ha 301 are indicated. V1647 Ori has an H-band magnitude of 11.96. 
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2 2.2 2.4 



Fig. 4. Ncar-IR spectroscopy of V1647 Ori taken on the NASA IRTF using SpcX on UT 22 February 
2007. Both the SXD and LXD grisms were used producing X-dispersed spectra of the 1 to 2.5 fxm and 2.5 
to 4.1 /im region, respectively. The top panel shows the complete spectrum while the middle and bottom 
panels show the J+H and K band spectral regions, respectively. The spectrum is calibrated in units of F^ 
such that the flux at 1.65 is 1.94 10"^'' W m~^ nm~^. The spectrum was normalized to unity at 2.2 /xm 
where the flux is 4.54 10~^^ W m~^ /im"^. 
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Fig. 5. — An expanded view of the 2.1-2.35 region of the near-IR spectroscopy of V1647 Ori. The spec- 

tmiri is calibrated in units of Fa and normalized to unity at 2.2 fim where the flux is 4.54 10""^^ W /im"-^. 
Note the presence of weak CO bandhead absorption at 2.294 /um. 
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Fig. 6. — An optical narrowband [S II] image of V1647 Ori, McNeil's Nebula, and the HH object HH 22 and 
HH 23. North is at the top, east is to the left. HH 23 is 155" north of V1647 Ori and appears to originate 
from it. A possible intermediate bow shock is marked as are the HH 22 jet, and the far-IR source detected 
by Spitzer, HH22-IRS (MuzeroUe et al. 2004). 
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Fig. 7. — A plot of the K-band spectrum of V1647 Ori (black) overlaid by model fits using the method of 
Prato, Greene, & Simon (2003). Three different model fits are shown for spectral types Ml V (orange), 
MO V (green), and K7 V (blue). These produce the best overall fit to the data. The derived Ay and r^ 
values for these models are 29 and 1.5, respectively. The poor fit longward of 2.3 /im is due to the presence 
of excess water vapor absorption from the circumstellar environment. 
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Fig. 8. — A plot of our adopted near-IR disk excess emission model. The red curve shows the shape of the 
near-IR disk excess emission of temperature 1000 K. In the blue curve, we have incorporated a rudimentary 
model for water vapor absorption by adding the shape of the water vapor features of an M7 V star. This 
was achieved hy i) removing all other spectral features from the M7 V data, ii) smoothing the data, and in) 
removing the M7 V continuum shape by subtracting an appropriate black-body curve. 
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Fig. 9. — The surface plot for model fits to the V1647 Ori 1-4 fim spectrum. In this plot we used an MO V 
star as representative of V1647 Ori to plot the derived Ay and Tdust values for different tk values ranging 
from 1.3 to 1.7. The surfaces for = 1-01, 2, and 3 are shown. The 3.0 surface corresponds approximately 
to la. The range of acceptable values of Ay and Tdust, AAy and ATdust, are ~2 magnitudes and ~250 K, 
respectively. 
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Fig. 10. — A similar plot to that shown in [5] but for a range of spectral types, K5 V to M2 V for fixed 
r/f=1.5. Again, the la surface corresponds to a range of AAy ^2 and ATdust '~250 K. 
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Fig. 11. — A plot of the observed (black) and best-fit model spectrum (red) using an MO V spectral type 
template star, Ay— 19, rx = 1.5, the water vapor model defined in the text (Fig. [5]), and incorporating a 
model ice absorption band from laboratory ices. The poor fit from 3.2 to 3.8 /im is due to the model not 
including the long- wavelength wing of the 3.1 ^m water ice absorption feature that has, as yet, an unknown 
cause. 
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800 1000 1200 1400 

Fig. 12. — A surface plot for the final adopted values for the best-fit model. Specifically, these are an 
MO V star, rj<-=1.5 and the water vapor model. The x^=10 surface corresponds approximately to 3a. 
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Fig. 13. — A near-IR J-H vs. H-K color-color diagram showing the published photometry for V1647 Ori 
from the pre-out burst and outburst phases together with our quiescent phase values. The photometric values 
are shown as solid circles and the ID numbers refer to the source and time of the observations. The legend 
shows the identification of the photometry, relating ID number to source of data and observing time. Point 
1 is from 2MASS data, point 2 from Reipurth & Aspin (2004) Gemini data, point 3 from the United States 
Naval Observatory telescope as presented in McGchcc ct al. (2004), points 4 and 6 arc from the Himalayan 
Chandra Telescope (HCT) presented in Ojha ct al. (2006), point 5 is from the Telescopio Carlos Sanchez 
(TCS) at the Teide Observatory presented by Acosta-Pulido ct al. (2007), and point 7 is from the data 
presented here. The colors of main sequence dwarfs (labeled V), giants (labeled III), and classical T Tauri 
stars (labeled CTTS) are also shown. The CTTS locus (the dashed line) is from Meyer ct al. (1997). 
Reddening vectors corresponding to Av=2G magnitudes are shown extending from the extremes of the main 
sequence/giant branch (solid lines) and the same reddening vector is shown extending from the extreme of 
the CTTS colors (dash-dotted line). 
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Fig. 14. — A plot of Av vs. Modified Julian Date (MJD) showing the change in Ay with time from pre- 
outburst to outburst, and into quiescence. The Ay values were derived from the J-H and H-K colors in 
Fig. [13] by dereddening onto the CTTS locus from Meyer et al. (1997). The dashed line represents the Ay 
value derived from the 2MASS J, H, and K' colors which we adopt as our pre-outburst value. The outburst 
values are those from February 2004 to February 2006 while the quiescent value is from our February 2007 
data. Clearly, in early 2006 V1647 Ori reverted to its pre-outburst color. The difference in Ay between the 
pre-outburst and outburst values is AAy ^^8.5 magnitudes. 
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Fig. 15. — A comparison of our derived stellar luminosity L*=5.2 L0 and effective temperature Te//=3800 K 
with evolutionary track models from Baraffe et al. (1998). The location of V1647 Ori is indicated by a filled 
star with associated uncertainties displayed. The location of V1647 Ori in this plot suggests it is a star of 
mass ~1 Mq with an age younger than <0.5 Myr. 
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Fig. 16. — A similar plot to that shown in Fig. [15] but for the evolutionary tracks models of Siess et al. 
(2000). Again, V1647 Ori is indicated by the filled star with associated uncertainties displayed. The location 
of V1647 Ori in this plot suggests a stellar mass of ~0.55 Mq and an age of <0.5 Myrs. 
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Fig. 17. — The spectral energy distribution (SED) of V1647 Ori. Plotted are the pre-outburst SED (open 
squares) and the outburst SED (filled dots) both from Andrews, Rothberg, & Simon (2004), and the quiescent 
phase SED from our data (open stars). Uncertainties on the quiescent phase data are indicated. Note that 
the pre-outburst and quiescent phase SEDs are very similar. 
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Fig. 18. — The 8 13 /iin spectrum of V1647 Ori taken with the Gemini MicheUe imaging spectrometer on 
UT April 5, 2007. The spectrum was ratioed with a telluric standard star and calibrated using the Michelle 
photometry shown in Table 1. Hence, the flux at 11.2 /um is 0.23 Jy. The dashed line is a linear interpolation 
between 8 and 13 /um (see text for details). 
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Fig. 19. — A plot of the 8-13 silicate absorption band optical depth extracted from the observed spectrum 
shown in Fig. [TH]by subtracting the observed the spectrum of V1647 Ori from the UT April 5, 2007 data. 
The silicate absorption flux at each wavelength was subtracted from the flux from the linear interpolation 
as in Quanz ct al. (2007). The resultant profile was then converted to optical depth, r. 
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